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We investigate cosmological effects concerning the late-time entropy production due to the decay 
of non-relativistic massive particles. The thermalization process of neutrinos after the entropy 
production is properly solved by using the Boltzmann equation. If a large entropy production takes 
place at late time t~ 1 sec, it is found that a large fraction of neutrinos cannot be thermalized. This 
fact loosens the tight constraint on the reheating temperature Tr from the big bang nucleosynthesis 
and Tr could be as low as 0.5 MeV. The influence on the large scale structure formation and cosmic 
microwave background anisotropics is also discussed. 
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It is usually believed that thermal radiation dominates 
the energy density of the early universe after the reheat- 
ing process of the primordial inflation. At least, the uni- 
verse is expected to be radiation-dominated before the 
big bang nucleosynthesis (BBN) epoch (t ~ 1 sec) oth- 
erwise the abundances of the synthesized light elements 
( 4 He, 3 He, D and 7 Li) do not agree with observations pi . 
However, it is uncertain that the universe is radiation 
dominated before BBN epoch. In fact, particle physics 
models beyond the standard one predict a number of new 
massive particles some of which have long lifetimes and 
may influence the standard BBN scenario. Since the en- 
ergy density of such massive non-relativistic particles de- 
creases more slowly than that of the radiation as the uni- 
verse expands, the universe becomes matter-dominated 
by those particles until they decay. When they decay into 
ordinary particles, the large entropy is produced and uni- 
verse becomes radiation-dominated again. We call this 
process "late-time entropy production" . 

One can find some interesting candidates for the late- 
time entropy production in models based on supersym- 
metry (SUSY). In local SUSY (i.e. supergravity ) theo- 
ries [|| there exist gravitino and Polonyi field || which 
have masses of - O(100GeV - lOTeV). Since grav- 
itino and Polonyi field interact with other particles only 
through gravity, they have long lifetimes. For example, 
the Polonyi field with mass ~ 10 TeV quickly dominates 
the energy density of the universe because this field can- 
not be diluted by usual inflation and decays at BBN 
epoch. It is also known that superstring theories have 
many light fields such as dilaton and moduli which have 
similar properties as the Polonyi field. 

When one considers the late-time entropy production, 
reheating temperature Tr is usually used as a param- 
eter to characterize it. The reheating temperature Tr 
is determined from T = 3H(Tr) where T is the de- 
cay rate and H(Tn) is the Hubble parameter at the de- 
cay epoch. Since the Hubble parameter is expressed as 
H = y/g^/MTl/Ma, where <?» is the number of mass- 



less degrees of freedom (= 43/4) and Mq is the reduced 
Planck mass (= 2.4 x 10 18 GeV), the reheating tempera- 
ture is estimated as T R = 0.554^f^Mc. 

As mentioned above, the stringent constraint on the 
late-time entropy production or reheating temperature 
comes from the consideration of BBN. The long-lived 
massive particles which are responsible for the late-time 
entropy production should decay early enough to make 
the universe to be dominated by thermal radiation be- 
fore the BBN epoch. To establish the thermal equilib- 
rium, the decay products should be quickly thermalized 
through scatterings, annihilations, pair creations and fur- 
ther decays. Almost all standard particles except neutri- 
nos are thermalized very soon when they are produced 
in the decay and subsequent thermalization processes. 
Neutrinos can be thermalized only through weak inter- 
action which usually decouples at a few MeV. Thus, the 
thermalization of neutrinos is most important to obtain 
constraints on the reheating temperature. However, the 
thermalization of neutrinos has not been well studied and 
people have used various constraints on the reheating 
temperature between lMeV and lOMeV. Therefore, in 
this letter, we will obtain the constraint on the reheat- 
ing temperature by using the neutrino spectrum obtained 
from numerical integration of a set of Boltzmann equa- 
tions together with full BBN network calculations. 

Another interesting constraint may come from observa- 
tions of anisotropies of the cosmic microwave background 
radiations (CMB). It is known that the CMB anisotropies 
are very sensitive to the time of matter-radiation equal- 
ity (see e.g., |Q) . When the reheating temperature is 
so low that sufficient neutrinos cannot be thermally pro- 
duced, the radiation (= photons + neutrinos) density 
becomes less than that in the standard case, which may 
give distinguishable signals in the CMB anisotropies. In 
this letter, we use the effective number of neutrino species 
N* as a parameter which represents the the energy den- 
sity of neutrinos defined by N° s = ^ p Ui /p s td, where 
i = v e ,v^,v T , and p s td is the neutrino energy density in 
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the standard case (i.e. no late-time entropy production). 

First, let us discuss the neutrino spectrum. When a 
massive particle <j> which is responsible for the late-time 
entropy production decays, all emitted particles except 
neutrinos are quickly thcrmalized and make a thermal 
plasma with temperature ~ Tr. If the reheating temper- 
ature is high enough (Tr ^> 10 MeV), there is no question 
about the neutrino thermalization. For relatively low re- 
heating temperature (Tr < lOMcV), however, neutrinos 
are slowly thermalized and may not be in time for the 
beginning of BBN. We assume that the decay branch- 
ing ratio into neutrinos is very small and that neutrinos 
are produced only through annihilations of electrons and 
positrons, i.e. e + + e~ — > z/, + Vi(i = e, /i, r). The evolu- 
tion of the distribution function fi of the neutrino i/j is 
described by the Boltzmann equation: 



dfi(pi,t) _ Hp d fi(Pi,t) 



at 



Cf + Cf 



(1) 



where pi is the momentum of Vi and Cf (Cf ) is the col- 
lision term for annihilation (scattering) processes. Here 
we consider the following processes: Vi + V{ «-» e + + e~ 
and Vi + <-> z/j + e ± . We do not include the neutrino 
self interactions which may not change the result in the 
present paper, since the neutrino number densities are 
much smaller than the electron number density at low 
reheating temperature. 

Here we have treated neutrinos as Majorana ones (i.e., 
v = v). Note that our results are the same for Dirac 
neutrinos. The collision terms are quite complicated and 
expressed by nine dimensional integrations over momen- 
tum space. However, if we neglect electron mass and 
assume that electrons obey the Boltzmann distribution 
e _p / T , the collision terms are simplified to one dimen- 
sional integration forms. Because the weak interaction 
rate is small at T < 0.5 MeV, neglecting the electron 
mass changes the result little. Then, Cf is given by || 



C? 



dp': 



yfiavUf^fM) - U q {Pi)feM)\ (2) 



where f e q(= l/(e Pi ^ T + 1)) is the equilibrium distribu- 
tion and (crv)i is the differential cross sections given 
by Me = (4G|/(9tt))((C v + l) 2 + [C A + l) 2 )pp> , 
and (av)^ T = (AG 2 f /(9tt))(C^ + C\)pp' . Here G F is 
the Fermi coupling constant, and Ca — — 1/2, Cy = 
-1/2 + 2 sin 2 (9 W - Weinberg angle). 

As for scattering processes, Cf is expressed as 



Cf 



2G ^-(C 2 +Ci)x 
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Here (C 2 + C\) is replaced by ((C v + l) 2 + (C A + l) 2 ) 
for i = e, and the functions F±, F 2 , Bi, B 2 are given by 

FripJ) = D{p,p>) + E{p,p')e-v' ' T , 
F 2 (p,p>) = D{p' lP )e^' T + E(p,p')e~P'/ T , 
B 1 ( P ,p')=F 2 (p'p), B 2 (p,p') = F 1 (p'p), (4) 

where 

D(p,p') = 2T i (p 2 +p' 2 + 2T(p-p') +4T 2 ), 
E(p,p') = -T 2 [p 2 p' 2 + 2pp'(p + p')T 

+2(p + p') 2 T 2 + 4(p + p')T 3 + 8T 4 ]. (5) 

Together with the above Boltzmann equations, we also 
solve the evolution of the densities of the massive particle 
4>, radiations p r (= p 7 + p e ±) an( i the scale factor a: 



dp,}, 
~dt 

dp r 



-Tpcf, - 3Hp<p, 
dp v 



(6) 
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where the neutrino density is given by p v 
Ei 1 /* 2 f d PiPiPifi(Pi)- 
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FIG. 1. (a) Evolution of the the energy density of z/ e (solid 
curve) and (dashed curve) for Tr = 2 MeV. (b) Distribu- 
tion of v e (solid curve) and (dashed curve) for Tr — 2 MeV. 
The dotted curve is the thermal equilibrium Fermi-Dirac dis- 
tribution. 

In Fig. |l| (a) we show the evolutions of p Va and p Ufi = 
p Vr for Tr = 2 MeV. Since the electron neutrinos interact 
with electrons or positrons through both charged and 
neutral currents, they are more effectively produced from 
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the thermal plasma than the other neutrinos which have 
only neutral current interactions. The final distribution 
functions f e and / M = f T are shown in Fig. |l| (b), from 
which one can see that the occupation numbers are close 
to equilibrium values at low momentum but they deviate 
significantly at higher momentum. 

In Fig. U we can see the change of N° s as a function 
of T R . ttT R >7 MeV, Nf is almost equal to three and 
neutrinos are thcrmalized very well. On the other hand, 
if Tr < 7 MeV, Nf becomes smaller than three. 

The deficit of the neutrino distribution influences the 
produced light element abundances. In particular, the 
abundance of the primordial 4 He is drastically changed. 
At the beginning of BBN (T - 1 MeV - 0.1 MeV) the 
competition between the Hubble expansion rate and the 
weak interaction rates determines the freeze-out value 
of neutron to proton ratio. After the freeze-out time, 
neutrons can change into protons only through the free 
decay with the life time r„. Since the left neutrons are 
almost included into 4 He, the primordial 4 He is sensitive 
to the freeze-out value of neutron to proton ratio. 




T R (MeV) 
FIG. 2. Nf as a function of T R . 

If the neutrino energy density gets smaller than that 
of the standard BBN (SBBN), Hubble parameter is also 
decreased. Then the (3 equilibrium between neutrons and 
protons continues for longer time, and less neutrons are 
left. Thus, the predicted 4 He is less than the prediction of 
SBBN. This effect is approximately estimated by: AY" ~ 
— 0.1(— Ap/p), where Y is the mass fraction of 4 He and 
p is the total energy density of the universe. 

In addition if the electron neutrino is not thermalized 
sufficiently and does not have the perfect Fermi-Dirac 
distribution, there are two interesting effects by which 
more 4 He are produced. The weak reaction rates are 
computed by the neutriono distributions which are ob- 
tained by solving Boltzmann equations. For example, a 
reaction rate at which neutron is changed into proton 
(ri + i> e — > p + e~) is represented by 



where Q — m n — m p — 1.29 MeV and K is a nor- 
malization factor which is determined by the neutron 
life time r n as K ~ (1.636r„) _1 . In this equation we 
can see that if the (anti) electron neutrino distribution 
functions decrease, the weak interaction rates also de- 
crease.0 First, when the weak interaction rate T n ^ p de- 
creases, the Hubble expansion becomes more rapid than 
that of the interaction rate earlier. Then the freeze-out 
value of neutron to proton ratio becomes larger than in 
SBBN and the predicted 4 He abundance becomes larger: 
AY ~ +0.15(-Ar„^p/r n „ p ). Second when the inter- 
action rates T n —> p at which neutrons are changed into 
protons become smaller, less neutrons can turn into pro- 
tons after the freeze-out time. Then the produced 4 He 
also becomes larger: AY" ~ +0.2(— Ar n _> p /r n _>£,). 

As for the observational abundances, we adopt the fol- 
lowing values. The primordial 4 He mass fraction Y is 
observed in the low metallicity extragalactic HII regions. 
Now we have two observational values, low 4 He and high 
4 He, which are reported by different groups. We take 
"Low 4 He" from Olive, Skillman and Steigman (1997) 
@, Y° bs = 0.234 ± (0.002) stat ± (0.005) syst . Recently 
Izotov et al. M claimed that the effect of the Hel stel- 
lar absorption which are not considered well in j|] is 
very important. We adopt their value as "High 4 He", 
Y obs = 0.244 ± (0.002) stat ± (0.005)^. 

The deuterium D/H is measured in the high redshift 
QSO absorption systems. Here we adopt the most reli- 
able data D/H = (3.39 ± 0.25) x 10~ 5 §. 

The 7 Li/H is observed in the Pop II old halo stars. 
We take the recent measurements ||] and adopt the ad- 
ditional larger systematic error, for fear there are under- 
estimates in the stellar depletion and the production by 
the cosmic ray spallation. Then we obtain: Logio( 7 Li/H) 
- -9.76 ± (0.012) stot ± (0m) syst ± (0.3) Qdd . 

In order to discuss how the theoretical predictions with 
the late-time entropy production agree with the above 
observational constraints, we perform Monte Carlo cal- 
culation and the maximum likelihood analysis [ for de- 
tails, e.g. see po| l. In Fig. |we plot the contours of the 
confidence level in the tj-Tr plane for (a) Low 4 He and 
(b) High 4 He. As we mentioned above, for Tj? > 7 MeV, 
the theoretical prediction is the same as SBBN. On the 
other hand as Tr decreases, Y gradually becomes smaller 
because the effective number of neutrino species Nf de- 
creases. On the other hand, for Tr < 2 MeV the effect 
that the weak interaction rates are weakened due to the 
lack of the neutrino distributions begins to be important 
and Y begins to increase as Tr decreases. For Tr < 1 
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*The weak interaction rates with (anti) electron neutrino 
in the final state slightly increase because the Pauli block- 
ing factor (1-fue) increases. However, both of the total weak 
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MeV, the weak interaction rates are still more weakened 
and Y steeply increases as Tr decreases because it is too 
late to produce enough electrons and positrons whose 
mass is about m e = 0.511 MeV. From Fig. | we can con- 
clude that T R < 0.5 MeV is excluded at 95 % C.L. Or 
we could rather say that Tr can be as low as 0.5MeV. 
Accordingly N^ s can be as small as 0.1. 
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FIG. 3. Contours of the confidence levels in (j), Tr) plane 
for (a) observational value of Low 4 He and (b) High 4 He. The 
inner (outer) curve is 68% (95%) C.L.. 

The late time entropy production is also constrained 
from the formation of the large scale structure of the 
universe. We can fit the large scale galaxy distribution 
if the "shape parameter" T s ~ 0.25 ± 0.05 Jn|. In case 
of the standard cold dark matter (CDM) scenario, T s = 
Qoh whose dependence is determined by the epoch of 
the matter-radiation equality. Here CIq is the density 
parameter and h is the non-dimensional Hubble constant 
normalized by lOOkm/s/Mpc. Since T s ~ 0.5 for so- 
called standard CDM model, it is known that this model 
is in serious trouble and one can achieve the desired fit 
with low density models where fig — 0.3. If the late time 
entropy production takes place, however, T s is modified 
as T s = 1.68fl Q h/{l + 0.23Nf). Therefore the constraint 
from the large scale galaxy distribution becomes much 
tighter with Nf < 3 (see Fig. 4). 
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FIG. 4. Contours of T s = 0.2 (bold), 0.3,0.4,0.5 and 0.6 
on the (fi 0l Nf) plane for h = 0.7. 

Finally we discuss the CMB constraint on Tr. With 
the present angular resolutions and sensitivities of COBE 
observation ]12|| or current balloon and ground base ex- 



periments, it is impossible to set a constraint on N£ s . 
However it is expected that future satellite experiments 
such as MAP @ and PLANCK @ will give us a useful 
information about . From Lopez et al.'s analysis [0, 
MAP and PLANCK have sensitivities that SNf > 0.1 
(MAP) and 0.03 (PLANCK) including polarization data, 
even if all cosmological parameters are determined simul- 
taneously (see also Fig. 5). From such future observations 
of anisotropies of CMB, it is expected that we can pre- 
cisely determine Tr. 
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and polarization (right top panel) of models with 7V° ff = 3, 2 
and 0.5. Bottom two panels show (Ce(N e s) — Ci(3))/Ce(3) 
with iV c c | = 2.9, 2.5 and 2 for CMB anisotropies (left bottom) 
and polarization (right bottom). 
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